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of a tri-nuclear molecular structure before the cold ternary fragmentation of 
252 Cf takes place, we solved the eigenvalue problem of a certain class of vibra- 
tions which are very likely to occur in these molecules. These oscillations are 
the result of the joined action of rotations of the heavier fragments and the 
transversal vibrations of the lighter spherical cluster with respect to the fission 

X' 

' axis. In the calculation of the interaction between the heavier fragments we 

took into account higher multipole deformations, including the hexadecupole 
one, and introduced a repulsive nuclear part to insure the creation of a po- 
tential pocket in which a few molecular states can be accommodated. The 
possibility to observe the de-excitation of such states is discussed in connection 
with the molecular life-time. 
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The scope of this letter is to extend our recent investigations on the molecular configura- 
tions in the binary cold fission and the ternary cold fission ||. Growing interest aroused 
in the last year due to the experimental indication of a long-living (> 10 _13 s) structure in 
the 10 Be accompanied ternary cold fission Q]. Very recently a molecular structure in which 
12 C plays the role of the light accompanying particle has been reported ||5[. 

In this letter we make the suggestion that the ternary cold fission of 252 Cf is a process 
consisting of two main stages : in the preformation stage a quasi-bound molecular structure 
is formed in which the heavier fragments are almost co-linear and the light particle (e.g. 
a, 3 Li , 10 Be and 12 C) which is responsible for the molecular bonding, is orbiting in the 
equatorial region. This is similar to the case encountered in Molecular Physics, where in a 
linear or nonlinear chain of three atoms, the central atom ensures two bondings with the 
eccentric atoms In the second stage the quasi-molecular state is decaying f7|. Our aim 
is to study the collective vibrations of the system before the decay takes place. 

Recently it has been advocated by us based on the concept of nuclear molecule ||, 
that for fragments emitted in the binary cold fission, with almost no excitation energy, a 
collective vibrational spectrum will show up as a consequence of small non-axial fluctuations 
at scission. Such a molecular spectrum can be achieved if the interplay between the Coulomb 
and the repulsive nuclear core on one hand and the attractive nuclear part on the other hand 
will produce a pocket in the interaction potential between the fragments IJ. In the case of 
di-nuclear systems it was shown that possible molecular collective modes can be associated to 
the elongation variable and rotational vibrations taking place perpendicularly to the fission 
T0|1 • The last type of modes, e.g. butterfly (bending) and anti-butterfly (wriggling) is 
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also believed to be responsible for the formation of angular momenta in fragments emerging 
in binary spontaneous fission P, H],P^|. 

In a previous paper the classical expression of the tri-nuclear Hamiltonian has been 
worked out for the case of the 96 Sr+ 10 Be+ 146 Ba molecule in terms of the Jacobi variables 
R, £ and the angular velocities u/ of the molecular frame . The equilibrium configuration 
was that of three aligned clusters, with the lighter in-between. In such a configuration the 



interaction between the heavy fragments is almost entirely given by the Coulomb term. 
However the interaction between the lighter fragment and the heavy fragments consists also 
of a noticeable nuclear component, which in fact is responsible for the nuclear bond. Like 
in the case of binary molecules, butterfly modes can occur, in which the fragments rotates 
in phase while the lighter fragment is approximately preserving its pole-pole configuration 
with the heavier fragments. 

The classical expression of the kinetic energy of the three-body system, after removing 
the center of mass contribution, is expressed as a sum of translational and rotational degrees 
of freedom : 

T = ^ 12 R 2 + ^(12) 3 £ 2 + -£U\J\U\ + -*W2.72W2 + -*u> 3 J 3 U 3 (1) 

The first term describes the relative motion of the di-nuclear sub-system (12) with reduced 
mass /ii2 = m\m2j{jn\ + 777.2), whereas the second one corresponds to the relative motion 
of the third cluster with respect to the heavier fragments center-of-mass with reduced mass 
A*(i2)3 — ( m i + rn 2 )m 3 / (m 1 + m 2 + 777,3). The vectors u>i t 2,3 denote the angular velocities of 
the rotational motion of the three clusters, referred to the laboratory frame, t u> being the 
transpose of uj. In this paper we consider a spherical light cluster and thence the last term 
in eq.(|l]) disappears. The inertia tensors J" i are defined in the intrinsic frame such that 
the only non-vanishing components are the first two diagonal terms, (3~i)ix = {$'1)22 = J%, 
the quantum rotation around the symmetry axis of any of the two heavier fragments being 
discarded. 

In what follows we are interested in studying the collective spectrum which develops 
upon constraining the tri-nuclear molecule to perform an oscilllatory motion similar to the 
valence angle bending in molecular physics and the butterfly (bending) modes in di-nuclear 
molecules, i.e. to perform small displacements from the equilibrium position which result in 
the decrease of the angle between the two valence bonds, $ = 7T — (pi — (p%, attached to the 
spherical light fragment 3 (ifi is the angle between the axis joining the two heavier fragments 
and the line joining the heavy fragment % with the light cluster). In the same time, since the 
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nuclear proximity forces have the tendency to keep constant the reciprocal distances and 
orientations of the heavy fragments with the light one, we exclude possible bond stretching 
vibrations. If the bond stretching is absent, then there will be a corresponding decrease 
in the distance between the heavy nuclei 1 and 2, when the bending angles (pi and p 2 are 
increasing. The quantitative translation of the above mentioned considerations provides us 
with a set of constraints between the variables of interest in this problem : i?-the distance 
between the centers of the two heavier fragments, £-the distance between the light cluster 
3 and the center-of-mass of the heavy fragments ensemble, and the small bending angles 
ifi, tp2- These last two variables are related between them, due to the assumption on the 
constancy of the pole-pole configuration between the light cluster 3 and the heavy fragments 
1 and 2 

V2 = WTW 3 E (2) 

where e — (pi. Consequently we obtain the following relations, which allows us to eliminate 
from the kinetic energy ([!]) the variables R and £ in favor of e 

R=(R 1 + R 2 + 2R 3 ) (l - If^e") (3) 



e=0 



where 



t A 1 (R 1 + R 3 ) - A 2 (R 2 + R 3 ) 

e ° = e(0) = A\ + A 2 (5) 



MM (R 1 + R 3 )(R 1 + R 2 + 2R 3 ) 2 

de 2 ) £=0 A, + A 2 (R 2 + i? 3 )(A(i?i + R 3 ) - A 2 (R 2 + R 3 )) 



(6) 



Note that the above expressions have been written in the second order in e. 

With the above choice the heavy fragments are constrained to rotate only around an axis 
perpendicularly to the axis joining their centers. This possibility is justified experimentally 
by the small forward anisotropy of the angular distribution of prompt 7 radiation. 

As well as the above approximations we consider that the nuclei, building-up the 
molecule, are not performing (3 or 7 vibrations. 



We define a molecular frame whose z-axis coincides with the fission axis and the three- 
body plane is choosen to coincide with the x — z molecular plane (in this way we eliminate 
the ^/-component of the Jacobi coordinate £). Also, some assumptions have been made for 
the Euler angles of the interacting deformed fragments (xi, <Pi, 4>i)- The angle has already 
been defined above : it describes the angle between the fragment % symmetry axis and the 
molecular z-axis and it is expressed in terms of e (see eq. @j ) . The geometry of our problem, 
with the heavier fragments symmetry axes lying in the same plane, makes the Euler angles 
X equal , i.e. xi — X2- They are combined in the variable 9 3 = (xi + X2V2 which measures 
the rotation of the tri-nuclear aggregate with respect to the fission axis. 

Following the standard procedures H8|JT3| the total kinetic energy (P can be expressed 
as a sum of three parts, the rotational energy T rot , the internal kinetic energy T int and 
the Coriolis coupling T cor . In terms of the time derivatives of the Euler angles, (#1, # 2 , #3), 
specifying the rotation of the molecular frame, the classical, rotational kinetic energy reads: 

Trot=\Y,9lfeji (7) 
Z ij 

where the only non-vanishing components of the rotational metric tensor g™ 1 are given by : 

9u = (Jo + Ji + J2) sin 2 9 2 - J 13 sin 29 2 cos 9 3 + (j t + J 2 f 1 + f 3 ") e sin 29 2 cos 29 3 

\ R 2 + R3 / 

#22 = ^0 + J\ + J2 

Ri + R 



9n ={jn- 2e cos 9 3 (j x + h ^ R * ) } sin #3 cos 9 2 
923 = [Jn - 2ecos^ 3 (Ji + J 2^^)} sin ^3 

9n = ~ { Jis cos 9 3 - e cos 29 3 (ji + J 2 ^^ ) } sin 9 2 (8) 

where J = [i\ 2 [R\ + R 2 + 2R 3 f + ^(i 2 ) 3 £o an d J13 = ^(12)3(^1 + ^3)^0 

The intrinsic kinetic energy will be comprised of "e" -vibrations and intrinsic rotations 
of the clusters: 

Tint = \ ( J i3 jRl ^ 3 +J1 + J2 ( ^^ ) ) - J i3 (sin 2 sin 9 A + cos #3^ 



There will be also a contribution from the Coriolis term 

T cor = [h + J 2 ^±3.) { £ ( s in 2 cos 20 3 a - sin 26 3 6 2 )6 3 + (sin 6 2 sin 26 z 6\ + cos20 3 2 )£r} 

(10) 

The next development of our considerations is facilitated by the peculiarities of the 
spontaneous fission of 252 Cf. Since the spin J of the mother nucleus is 0, the total helicity 
K, i.e. the projection of the total angular momentum on the fission axis is zero for both 
binary and ternary fragmentations and therefore 83 = 0. To simplify even further one 
choose the molecular frame such that #3 = 0. After quantizing the kinetic energy in three 
coordinates (e, 61,62) and neglecting terms multiplied by the non-diagonal matrix-element 
Ji + ^2 ^^3 ~~ -^13) which prove to be small in the resulting metric tensor, we arrive to a 
form of the kinetic energy in which the rotations are decoupled from the butterfly vibrations 

f _ h 2 f i d 2 8 

2( Jo + Ji + J2) U n 62 d6 2 + COt 2 d6x + dd 2 ) 2J £ de 2 ( } 

where 

Next we turn our attention to the computation of the potential. The total interaction 
energy is given by the sum 

y= E Vu(Ru) (13) 

The interaction between two clusters composing the giant molecule can be calculated as the 
double folding integral of ground state one-body densities pi(2)(f) of heavy ions: 

V(R) = Jdr.J dr 2 pi(r 1 )p 2 (r 2 )i;(a) (14) 



We employ the M3Y iViV effective interaction for the nuclear part of v as described in [|14] to 
which we add a repulsive core in order to take into account two major factors - the density 
dependence of the NN interaction and the Pauli principle, which are important at distances 
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corresponding to the overlap of the nuclear volumes. This choice is particularly useful for a 
molecular model in which the repulsive core prevent the re-absorption of the lighter fragment 
by the heavier one. 

We consider that the nuclei composing the giant molecule are in their ground state with 
known quadrupole fl 2 , octupole (3 3 and hexadecupole deformations /3 4 . 

From what has been said above the interaction between the light cluster 3 and the heavy 
fragments remains unmodified when the bonding angle is decreased. The reciprocal distances 
and orientations between the light cluster and the heavier fragments being freezed. Therefore 
its contribution to the total Hamiltonian adds only a constant term. On the contrary, the 
interaction of the two heavy fragments 1 and 2 depends on the butterfly angle e as can be 
seen from the multipolar expansion of the potential 



Vi2= E 7== — = igWUo)y H U-^,o (is) 

Ai,a 2 ,a 3 ,/i y(2Ai + 1)(2A 2 + 1) \R2 + Rs J 

For small non- axial fluctuations(bendings), the potential in the neighborhood of the scission, 
or "molecular equilibrium" point R = Ri + R 2 + 2R 3 , gets a simplified form, provided we 
keep terms up to the second power in angle : 

V 12 = V(Ro) + l -C £ e 2 (16) 



where the stiffness parameter reads 



Ai A2A3 



MAi + + ( A 2 (A 2 + + A 3 (A 3 + F A Rl + ^ 3 



R 2 + R3 V R 2 + R3 ) R 2 + R3 



Vx^AaAaC-Ro) 



Hence the quantized vibrational Hamiltonian of the giant tri-nuclear molecule acquires the 
form 

h 2 d 2 1 

H v jh = — - _ H — C F £ (18) 

2J £ de 2 2 v ; 

The spectrum of the butterfly vibrations is given then simply by 
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E £ =(n £ + ~^hu £ (19) 

where oj £ — \/%. 

V 

As one can see from Table I, in the tri-nuclear case the values of Hu £ are only slightly 
smaller if we take into account Nuclear + Coulomb forces compared to considering only 
Coulomb forces. This was to be expected since in this case the distance between the two 
heavier fragments is already some fm beyond the top of the barrier and therefore the Coulomb 
forces are clearly dominating the interaction. This has to be contrasted to the case when 
we remove the lighter cluster, i.e. we insert R3 = and A 3 =0 in (^^), and the two heavier 
fragments are brought in touch. In this case the nuclear contribution to the interaction 
significantly decreases the value of hu e . Notice also the variation with the type of nuclear 
collective flow reflected in the type of inertia moment, the smallest values for hio £ being 
obtained for rigid rotations. 

To have an ideea of the order of magnitude of tri-nuclear molecules life-times we evaluate 
the half-life of the corresponding di-nuclear configuration, when the lighter cluster is absent. 
Using the one-dimensional WKB formula 

rR 3 t 



A = — - exp { —2 



ca P \ ^ \ fc 2 

27T I jR t2 ^ fl 

for the decay rate of a metastable state of energy 

fuo R 1 



2/212 ' D - — ^ I \ (20) 



1 (dV 12 (R) 



! _ (21) 



2 2 \| /i 13 V dR 2 

in the potential pocket of depth D = VuiRmax) — Vu^Rmin) of the two heavier nuclei, we 
obtain the half-life: 

ln2 . . 

Ti/2 = — (22) 

The computed life-times for di-nuclear molecules have large values(^ 10~ 13 s, according 
to Table I) and since the light particle has the tendency to delay the penetration of the 
mutual barrier of the heavier fragments, due to the attractive interaction with the heavier 
fragments, it is then justified to expect even larger values of Ti/ 2 for tri-nuclear molecules. 
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We choosed in this paper a linear configuration for the giant tri-nuclear molecule. Ac- 
cording to very recent calculations of penetrabilities in alpha ternary cold fission this con- 
figuration is reached at a certain step of the tunneling process J/] . It could be possible that 
before the beginning of mutual penetration of the multi-dimensional barrier, the giant tri- 
nuclear molecule is found in a tri-angular quasi-equilibrium configuration. In a forthcoming 
paper we will present the calculations of the molecular collective spectra also for such a 
configuration but the general trends should be same : considering an adiabatic scenario for 
the decay, the collective molecular modes will adjust slowly to the value of the elongation 
variable. As we mentioned earlier, the life-times of such nuclear molecules are suspected to 
be larger than 10 _13 s, whereas the life-times of the first molecular states computed in this 
paper are expected to be r pa 1/uj £ > 10~ 22 s. The gammas coming from the de-excitation 
of these states should be observed before tunneling ended and therefore they should not be 
pronouncedly Doppler-shifted. 
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TABLES 

TABLE I. The quantum energy hoj £ of the butterfly mode in the case of tri-nuclear ( 10 Be and 
12 C accompanied ternary fission) and di- nuclear configurations. N + C signify calculations with 
nuclear and Coulomb forces, whereas C only with Coulomb forces. Three types of the inertia 
moments are considered: a) experimental , b) Irrotational Fluid and c) Rigid Rotator. On the last 
column we listed the half-life of the di-nuclear molecule. 



Splitting 



Inertia moment 

J = Jexp 
J = Jjp 

J = Jrr 



Tri-nuclear 
fiwf +c (KeV) ftcjf(KeV) 
1900.4 2040.4 
1691.2 1815.9 
1014.4 1089.1 



96 Sr+ io Be+ i46 Ba 

hto^ +c (KeV) 
1100.7 
880.9 
440.4 



Di-nuclear 
huf (KeV) 
3377.3 
2702.8 
1351.4 



T l/2(s) 



< 5.5xl0" 9 



Splitting 
Inertia moment 

J = Jexp 
J = J if 

J = Jrr 



96 



Sr+ 12 C+ 144 Xe 



Tri-nuclear 
tiLU? +c {KeV) ^f(KeV) 
1799.7 1917.3 
1611.5 1716.8 
996.0 1061.1 



flUJ 



N+C 



(KeV) 



791.4 
625.5 
312.75 



Di-nuclear 
fiwf (KeV) 
3459.6 
2734.5 
1367.25 



T l/2(s) 



< 9.6xl0~ 9 
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